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Abstract 

We study the thermal phase diagram of pure SU(3) lattice gauge theory with fun- 
damental and adjoint couplings. We improve previous estimates of the position of 
the bulk transition line and determine the thermal deconfinement transition lines for 
Nt = 2, 4, 6, and 8. The endpoint of the bulk transition line {(3f, (3a) = (4.00(7), 2.06(8)) 
improves upon earlier estimates obtained using smaller lattice sizes. For A^^^ > 4 the 
deconfinement transition line splits cleanly away from the bulk transition line. With 
increasing Nt the thermal deconfinement transition lines shift to increasingly weaker 
coupling, joining onto the bulk transition line at increasingly larger f3a in a pattern 
consistent with the usual universality picture of lattice gauge theories. We also discuss 
the possible consequences of an induced adjoint term from the fermionic determinant 
and determine this induced term numerically with a microcanonical demon method for 
two flavors of dynamical Wilson fermions. 



1 Introduction 



The phase diagram of fundamental-adjoint pure gauge systems, shown in Fig. |1], is 
considerably more complicated than the one with only a fundamental coupling. Early 
studies found that it has first order (bulk) transitions in the region of small /?/ |T|, . For 
SU(3) the purely adjoint system has a transition at Pa = 6.5(3). For /3a — > oo the system 
reduces to a gauge theory with a transition at /?/ = 0.67. In the Pa,Pf plane these 
transition points are extended to lines that merge into a single terminating line segment 
as shown. As discussed in Sec. ^ we place the endpoint at (/?/, Pa) = (4.00(7), 2.06(8)), 
considerably different from old results obtained on small lattices 0. We also find a 
roughly straight line of bulk crossovers extending beyond the endpoint. 

The non-trivial phase structure in the fundamental-adjoint plane, and in particular 
the critical endpoint of the transition line pointing toward the fundamental axis, has 
been shown to be associated with, or even responsible, for the dip in the discrete /?- 
function of the theory with standard Wilson (fundamental) action, which occurs in the 
region where the bulk crossover line crosses the fundamental axis. The bulk transition 
might also mask the thermal deconfinement transition, at least for lattices with small 
temporal extent Nt. Indeed, from early simulations by Batrouni and Svetitsky and 
by Gocksch and Okawa IQ we know that for pure gauge SU(4) theory, for which the 
bulk transition crosses the fundamental axis, the Nt = 2 temperature transition joins 
the bulk line, but that for A^^^ = 4 the bulk line and the thermal line exist separately. 

In a recent paper Gavai, Grady and Mathur showed that at nonzero temperature the 
deconfinement transition in pure gauge SU(2) survives for positive adjoint couplings 
and is connected to the end point of the bulk transition line ||^. They were not able to 
distinguish the first order bulk line from the Nt = 4 and 6 thermal transition lines in 
their simulations. Universality in the continuum limit requires that as A^^^ is increased, 
the thermal phase transition line shifts to weaker coupling, so that approaching the 
zero-temperature weak-coupling limit along any line in the fundamental-adjoint cou- 
pling plane leads to the same low temperature, confined theory. Thus the thermal 
phase transition lines could not remain anchored to a bulk transition line for higher 
Nf. More recently Mathur and Gavai [^] raised doubts about the very existence of a 
bulk phase transition in SU(2), characterizing even the first-order portion of the line 
as a thermal deconfinement transition that is displaced toward weak coupling with 
increasing Nt. Clarification is obviously needed. 

Puzzled by the finding of Ref . , but also motivated by some unexpected results in 
the study of the high temperature behavior of full QCD with Wilson fermions that 
might be explained by the fact that the fermion determinant induces, among others, 
an effective adjoint coupling, we decided to study pure SU(3) gauge theory with a 
fundamental-adjoint action at finite temperature in more detail. We show that the 
first order bulk transition in pure gauge SU(3) in the fundamental-adjoint coupling 
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Figure 1: The phase diagram for the fundamental-adjoint system. 

plane indeed separates from the thermal deconfinement transition provided Nt is made 
large enough. 

According to universality arguments QCD with only two massless quarks has a 
second order thermal deconfining phase transition, which becomes a smooth crossover 
if the quarks are massive. Lattice QCD simulations with two light flavors of staggered 
quarks support this scenario. However, with two flavors of Wilson quarks recent sim- 
ulations indicate that at large values of the hopping parameter the phase diagram is 
more complicated 0: for A^^ = 6 and at hopping parameter k ~ 0.19 the system has 
a strong first order-like transition, where the plaquette expectation value has a sharp 
discontinuous jump. However, the average Polyakov loop remains small at this point, 
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and starts to increase only when (3 is considerably larger. This behavior is strongly 
reminiscent of separated bulk and thermal transitions, as in the pure gauge SU(3) 
model with mixed fundamental-adjoint action. In fact, these systems resemble each 
other so much that one is tempted to assume that the dynamical reason for these 
transitions is the same: the dynamical Wilson fermions induce a strong adjoint term 
in the effective pure gauge action. We therefore set out to measure such an induced 
adjoint term employing a microcanonical demon method. 

2 The fundamental— adjoint action 

Let us consider the action S for SU(A^) lattice gauge theory with fundamental and 
adjoint couphng: 

S = /3/E[l - ^ReTr^f^p] - ]^Tr^4Tr/t/p]. (1) 

Here Tr/- = Tr is the trace in the fundamental representation. Up is the path ordered 
product of the link matrices along the elementary plaquette P. Using Tr^f/ = |Tr/f/p — 
1 for the trace in the adjoint representation, the parametrization (|1|) is connected to 
the alternative one, used for general one-plaquette actions 0, 

^ = /^/ E[l - ^ReTr;[/p] + - ]^^Tr„f^p] (2) 

by 

Pa = J^^f^'a- (3) 

The A-parameters of the standard Wilson and the fundamental-adjoint action are 
connected by 

with bo the one-loop coefficient of the /5-function. Therefore, the equivalent Wilson 
coupling, in perturbation theory, is 

/3».^/3, + 2A-^^, (5) 

and lines of constant Pw in (|^) should represent fines of constant physics (up to two 
loop and non-perturbative corrections). Of course, at the relatively large couplings 
that we will be working with, the relation (^) is not very reliable. We shall try to 
improve the relation by using tadpole improved perturbation theory |TD|, where a 
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factor Mo, with Mq = Trllp/N, is taken out of every link field U^. In the relation (H) 
tadpole improvement amounts to the replacement Pa — utPa with uq computed in 
the equivalent Wilson theory. Compensating for the Uq as computed in perturbation 
theory then gives, instead 

Pw=Pf + 2PaUt - • (6) 

3 The phase diagram in the fundamental— adjoint 
coupling plane 

Monte Carlo simulations with the mixed action described above were carried out using 
a standard Metropolis updating scheme. First we found it necessary to repeat the 
early calculation of Bhanot to locate the endpoint of the bulk phase transition more 
accurately. This was done in a series of simulations with 8^ lattices by measuring the 
discontinuity in the average fundamental plaquette Ef = (TiUp) and adjoint plaquette 
Ea = (|Tr[/pp/3) along the first order line. Results are shown in Table |l|. To support 
the identification as a bulk transition Fig. ^ shows Ef, obtained from various lattice 
sizes, as function of (3f for (3a = 2.25 We note that Ef jumps for all lattices at the same 
point. 



(3a 


Pf 


AEf 


AEa 


3.0 


3.27 


0.656( 2) 


0.464(1) 


2.75 


3.47 


0.547( 2) 


0.398(2) 


2.5 


3.66 


0.424( 5) 


0.315(4) 


2.35 


3.77 


0.336( 9) 


0.250(7) 


2.25 


3.85 


0.252(10) 


0.193(8) 



Table 1: The jumps in in the average plaquette and average adjoint plaquette values 
across the bulk transition line. 



The location of the phase transitions for the different values of (3a considered was 
first obtained with runs from mixed starts. At the two (3f values that straddle the 
transition we then made runs with hot and cold starts to obtain the gaps in Table |l|. 
When we found metastable states in both cases we took the larger gap measured. In 
any case, the errors listed in Table |l] are statistical only. We estimate systematic errors 
to be of order 0.010 to 0.020 in AEf. 
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Figure 2: Ef as function of /?/ for various lattice sizes, as indicated, at = 2.25. 

The endpoint of the bulk transition line is then estimated by fitting to 

AEf = c{Pa- Plf. (7) 

Such a fit works very well, as shown in Fig. ^, having = 1-16 for 2 degrees of freedom. 
In this region of the parameter space the bulk transition line is essentially straight, and 
from the endpoint value (3* = 2.06(8) we infer the location of the endpoint of the bulk 
transition line as 

= (4.00(7), 2.06(8)), (8) 

shown with our improved location of the bulk phase transition line in Fig. |l]. The 
error of 0.08 in /?* includes our estimate of the systematic uncertainties in the gaps. 
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resulting as mentioned above from not obtaining them exactly at the critical coupling. 
The statistical error from the fit to the data in Table |l]is only 0.03. The best fit values 
for the other parameters are c = 0.68(1) and p = 0.58(3). 
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Figure 3: The average plaquette gap AEj as function of Pa together with the fit to 
determine the endpoint P*. 

4 Results at nonzero temperature 

As in pure SU(2), the SU(3) pure gauge deconfinement transition extends to non- 
zero Pa values, since the global Z(3) symmetry is not affected by the addition of an 
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adjoint coupling. As we increase the adjoint coupling the transitions for different Nt 
move closer together, and closer to the so called crossover region from strong to weak 
coupling behavior, which can be regarded as an extension of the bulk phase transition 
line. 

To explore the high temperature behavior of the mixed action theory, we studied the 
deconfinement transition on lattices with temporal extent Nt = 2, A, 6 and 8. Note 
that in the first 3 cases we limited ourselves to lattices with aspect ratio Ns/Nt = 2. 
That proved sufficient to find the deconfinement transitions with sufficient accuracy. 
The thermal deconfinement transition point is defined as the location of maximum 
slope in the modulus of the Polyakov loop (|-P|). All transition points we obtained are 
shown in Fig. ^ and Table |^. For Nf = 2 the = deconfinement transition continues 
smoothly into the bulk transition, which is shifted significantly from its large Nt value, 
as Pa increases. For Nt = 4 at f3a = 2.0 the deconfinement phase transition joins and 
is indistinguishable from from the bulk transition, just as was found for SU(2) [Q. But 
for Nt = 6 at Pa = 2.0, there is a clear separation, with the bulk crossover occurring at 
Pf = 4.035(5) and the thermal deconfinement transition at Pf = 4.07(2). At Pa = 2.25, 
even for the 12^ x 6 lattice, bulk and deconfinement transition coincide, both occurring 
at pf = 3.850(5). 



Pa 


Pf 


Nt = 2 


Nt = A 


Nt = Q 


Nt = 8 


0.0 


5.0941(4) 


5.6925(2) 


5.8941(5) 


6.001(25) 


0.5 


4.75(5) 


5.25(5) 


5.425(20) 




1.0 


4.4(1) 


4.85(5) 


4.96(2) 




1.5 


4.1(1) 


4.45(5) 


4.525(20) 


4.58(2) 


2.0 


3.75(5) 


4.035(5) 


4.07(2) 


4.135(10) 


2.25 




3.845(5) 


3.850(5) 


3.89(1) 


2.5 


3.45(5) 






3.660(5) 



Table 2: The deconfinement transition points. 

At larger Nt the small size of the order parameter required additional care in locating 
the phase transition. We carried out a finite size analysis with simulations on N^ x 8 
lattices with Ns = 8,12,16. In the confined phase {\P\) should vanish as l/\/Vs with 
increasing spatial volume, while in the deconfined phase it should extrapolate to a 
nonzero value. To treat both cases identical we made fits to {\P\) = Pod + cj \fVs- The 
finite size analysis is shown for Pa = 2.25 in Fig. ^. Plots for other Pa are very similar. In 
this way we place the Nf = 8 deconfinement transition for Pa = 2.0 at Pf = 4.135(15), 
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Figure 4: The phase diagram together with the thermal deconfinement transition points 
for Nt = 2, 4, 6 and 8. The lower plot shows an enlargement of the region around the 
end point of the bulk transition. 



clearly separated from the bulk crossover as well as from the Nf — 6 deconfinement 
transition. For = 2.25 the deconfinement transition occurs at Pf — 3.89(1), still 
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separated from the bulk transition at /3/ = 3.850(5), whereas at = 2.50, no clear 
separation is visible. 
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Figure 5: {\P\) at = 2.25 for Nf = 8 and A^^ = 8, 12 and 16 and the extrapolation 
to Ng — oo. 

We conclude that with increasing Nt the deconfinement transition line is displaced 
toward weaker coupling, joining onto the bulk transition at larger and larger values 
of Pa, consistent with universality. However, trying to see low temperature continuum 
physics at larger values of the adjoint coupling requires larger lattices to avoid strong 
violations of asymptotic scaling associated with the bulk transition line. 

Having found non-perturbatively lines of constant physics, here defined as lines of the 
thermal deconfinement transition for fixed N^, we can compare them to the perturbative 
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predictions obtained from eq. Not surprisingly, the perturbative prediction is poor. 
Using, instead, the tadpole improved prediction eq. changes the perturbative result 
in the right direction, but still fails to give the correct lines of constant physics, as 
shown in Fig. |^. It is interesting to notice, though, that Ef = (TiUp) is almost 
constant along the deconfinement transition lines. Therefore, equal effective coupling 
(3eS = 2/(1 — -Ej/S) describes lines of constant physics remarkably well. 




3 4 5 6 

Figure 6: Lines of constant physics as predicted by perturbation theory (dotted lines) 
and tadpole improved perturbation theory (dashed lines) together with the deconfine- 
ment transitions for Nt = 2, 4, 6, and 8. 
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5 Couplings induced by the fermionic determinant 



Full QCD simulations with two flavors of staggered fermions support the view that 
there is no thermal phase transition at positive quark masses, and a second order 
transition at zero quark mass |TT|] (a scenario also supported by universality arguments 



0). However, the phase diagram with dynamical Wilson fermion has turned out to be 
complicated for large values of the hopping parameter. It was found that for Nt = 6 
and hopping parameter k ~ 0.19, there is an apparent first order transition leading to 
a jump in the average plaquette. At slightly higher Pf the Polyakov loop expectation 
value (ReP) remains small in a manner remarkably similar to our findings for the region 
separating the bulk and thermal phase transition in the mixed action pure gauge theory. 
It is tempting to speculate that the first order phase transition seen with Wilson quarks 
is a bulk transition related to the first order pure gauge adjoint bulk transition. 

It has been shown, that the location of the deconfinement transition with staggered 
fermions can be surprisingly well explained by a change in the fundamental coupling 
induced by heavy fermions Let us see how new terms in the gauge action are 

induced by Wilson fermions. The discretized action 5*0 for n f flavors of Wilson fermions 
is ^ 

So = S, + j:Y.i'iMnm[UWm, (9) 
f=l n,m 

where 

) (10) 

and Sg is the pure gauge action. The fermions can be integrated out and the effective 
gauge action S^s can be obtained by a hopping parameter (k) expansion. 

^eff = Sg-UfTilnMp] 

= Sg + nfY. '^'^'"yj^Tr (jlil ± 7,)) • {TtUc + Trt/t ) , (11) 

where the sum is over all (unoriented) closed loops C and l[C] is the length of the loop. 
For small hopping parameters k, the shift for the fundamental coupling comes from 
loops around a single plaquette and is given by A/5/ = n/l^"'. The shift in the adjoint 
coupling comes from closed loops winding three times around a plaquette. For small 
K it is Aj3a oc K^^. Similarly the staggered fermion action induces an adjoint coupling. 



Indeed with eight flavors, a bulk transition has been found [13[, but studies have not 
yet been done to determine whether there is evidence for a similar separation of the 
bulk transition and the thermal crossover. 

The hopping parameter expansion outlined above is accurate only when k is very 
small. Results with much broader validity range can be obtained with the heavy quark 
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perturbation method. This method was used by Hasenfratz and DeGrand ||T2| to 
calculate A/3/. They found shifts A/3/ that agreed very well with MC data. However, 
such a heavy quark perturbation calculation has not been performed for A/3a. 



5.1 Demon Algorithm 



We use the micro canonical demon method p^ , |15|, |1^ to project out the induced 
effective gauge action from full fermionic simulations. Our ansatz for the effective 
action is given by eq. (||), with a priori unknown coupling constants Pi = /?/ or 
Pa- For both of the coupling constants we introduce a demon, which is a real- valued 
action variable with < Di < Df^^^. Effectively, the fermionic action is used as a 
heat bath to thermalize the demons: starting from a configuration generated with 
the original action, we update the system microcanonically with demons. During the 
microcanonical update Si + Di = const, separately for both fundamental and adjoint 
parts of the action. After the microcanonical update, the old gauge configuration is 
discarded, and the demon update is started again from a new gauge configuration, 
while preserving the demon values. The principle of the method is shown graphically 
in fig. 0. 



update with Sq 
O O 



copy 



O 



— " discard 

microcanonical 
update 

copy 



O 



O gauge configurations 
• demons 



Figure 7: Flow diagram of the demon algorithm. The configurations are generated 
with the original fermionic action 5*0 (top) eq. and copied for the microcanonical 
demon update. At the end of the microcanonical update, the demon values are copied 
for the next update phase (bottom). 

The microcanonical update phase is performed with a Metropolis algorithm: the 
proposed new gauge link is accepted only if the demons can give or take the amount 
of energy needed: 

AS^ > Di and - ASi < - A (12) 

for both i = f or a. Because both of the demons have to accept the update step for 
it to be performed, the demon distributions are correlated. After the whole update 
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procedure is repeated many times, the demons attain an equilibrium distribution and 
we can measure (-Dj). The induced effective couphngs can be solved from equations 



^^'^ ^ A " exmDf-) - 1 ■ ^^^^ 

It is not necessary to limit the demon action from above if the expected value of the 
corresponding coupling constant is large. In our case, we had to limit only the adjoint 
demon values from above; for the fundamental demon Z}™'^^ = oo and in this case 
eq. (|T3|) reduces to j3f = l/{Df). 

One drawback of the demon method is that the final results can depend on the 
details of the update procedure. Some of the factors affecting the results are the 
range and shape of the distribution from which the new gauge matrix is chosen in the 
Metropolis update (true even if the detailed balance is always satisfied), the acceptance 
rate, and the 'length' of the demon update phase: one can update each configuration 
microcanonically until the demons and the systems are properly thermalized, or one 
can stop the update after only one - or even partial - update sweep. These effects 
are easy to observe in some simple test models where the energy distributions can be 
exactly calculated. This does not mean that some methods are correct, some incorrect; 
different methods only correspond to different projections from the original action to 
the functional space spanned by the ansatz for the effective action. Nevertheless, one 
should not ignore this problem when using the demon method. Furthermore, let us 
note that if the effective action is completely equivalent to the original action, the above 
effects vanish and the demon method yields a unique Conversely, if the effective 
action is close to the true action - which should be the case for all reasonable effective 
actions - one can expect that the differences between the methods will be small. 

The autocorrelations between configurations can introduce further systematic errors 



T6|] . These errors are unrelated to the errors above. However, in most cases the 
differences are expected to behave as 1/V as ^ oo (the heat capacity of the system 
:§> the heat capacity of the demons). 

In our case we have a dramatically truncated effective action, and a priori we do 
not have any guarantee of the 'goodness' of the ansatz. The results were checked by 
using different update schemes: we performed 1 or 20 trajectories with the fermionic 
action between the microcanonical update phases, and 1/8, 1 or 50 microcanonical 
update sweeps for each gauge configuration. In all our tests the possible differences 
were completely overwhelmed by the statistical noise, indicating that the method is 
quite robust to the accuracy we reached. The results given below are all calculated by 
performing one microcanonical update sweep for each gauge configuration and after 
each fermionic trajectory. 



13 



O Q □ O 

5.5 ....dm. °^ 



o mo 

o □ 



5.0 







f 



I 



4.5 



cJ2 ^ 



4.0 



• • • 



0.10 0.12 0.14 0.16 0.18 0.20 0.22 

K 



Figure 8: The location of the dynamical Wilson fermion runs in the {f3'j, K)-plane 
(filled circles). The open circles indicate the location of the A^^ = 4 finite tempera- 
ture crossover, the squares the Nt = 6 phase transition, and the diamonds the zero- 
temperature Kc- The data is from ref. [j^ and from references therein. 

5.2 Simulations and Results 

We tested the demon method by applying it to pure gauge fundamental-adjoint simu- 
lations. On a 4'^-lattice we simulated the system at couplings (/3/,/3a)o = (3.6, 1.8) and 
(4.0,2.0). The induced couplings were (3.594(6), 1.812(17)) and (4.017(27), 1.92(6)), 
respectively, compatible with the input values. The latter coupling pair is very close 
to the bulk transition line in the (/5/, /5a)-plane. When we used an effective action 
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Figure 9: The induced Small crosses on horizontal dashed lines indicate the sim- 
ulation {Pj, k), open circles and black dots the measured /3/-values using the effective 
fundamental-adjoint and fundamental only gauge action, respectively. The length of 
the vertical bars gives the magnitude of the induced A/5/. 

consisting only of the Sj-part, /?/ was 4.688(9) and 6.044(12); the latter value is very 
close to the extension of the bulk transition line to the = -axis. 

If the Wilson fermion action induces a strong adjoint coupling giving rise to a bulk 
fundamental-adjoint transition, one should be able to observe the induced coupling 
already in small volumes. We performed simulations on 4^^ lattices with 28 different 
{(3p k) pairs. Fig. |^ shows the location of all runs performed. In Figs. § and we 
show the measured /5/ and /5a calculated with /3j = 4.0, 4.9, 5.46, and several k values. 
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Figure 10: The induced The scale on the right gives the magnitude of the induced 
I3ai measured as the vertical distance between the plot symbols and dashed horizontal 
lines. Points below the dashed lines indicate negative (3a- 



and in Fig. |TT| the results with constant n = 0.19. In Figs. and IT^ we also show jSf 
when Ses = PfSf only. 

When K is small, the quarks are very massive and the induced couplings are quite 
small (left side of Figs. ^ and 0). When k is increased, we approach the critical line 
where mq — >■ 0, and the fermionic contribution to the action becomes more significant. 
This is clearly visible as an increase in /?j in Figs. ^ and[Tl|. The critical values of n are 
approximately 0.16 {Pj = 5.46), 0.19 (4.9) and 0.22 (4.0). We observe no significant 
increase in Pa] on the contrary, when Hc is approached, the induced /?„ becomes slightly 
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Figure 11: The induced Pf and Pa, when k 



0.19. 



negative! The minor role of the adjoint action is also evident from the fact that /?/ 
remains practically the same whether we use the S'^-term of the effective action or not. 
We also checked the results with a few simulations on 6^ lattices with similar results. 
In Fig. |12|we compare A/5/ to the predictions of ref. [|12|, as a function of quark mass 



. The agreement is very good, especially when = 4.9. 
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Figure 12: Comparison of the analytical result for A/5/ |jl2|] by Hasenfratz and DeGrand 
(solid line) to the MC data. 

6 Conclusions 

We have shown that the fundamental-adjoint pure SU(3) gauge theory behaves as 
expected near the first order bulk transition. The thermal deconfinement transition 
lines join onto the bulk transition lines for larger and larger For Nt = 2 the thermal 
transition continues smoothly into the bulk transition which is, however, shifted from 
its location for larger Nf. The Nt = 4 thermal transition line joins the bulk transition 
line very close to the end point at Pa ~ 2.0, for Nt = 6 at Pa ^ 2.25 and for A^j = 8 at 
Pa ~ 2.5. Before joining the bulk transition line, the thermal transition line for a larger 
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Nt is to the right (at larger (3f) than for a smaller Nt. This finding is in agreement 
with what is expected from the usual universality picture of lattice gauge theories. 

We speculated that the bulk transition observed with Wilson fermions for Nt = 6 
and large k that looked very reminiscent of the situation seen in the pure gauge 
fundamental-adjoint action model, might be caused by an induced adjoint coupling. 
However, in measurements with the microcanonical demon method we did not observe 
that Wilson fermions induce any significant adjoint term in the pure gauge effective 
action, while the induced fundamental term is very well described by analytical calcu- 
lations. It is thus very improbable that the transition observed in iVt = 6 Wilson ther- 
modynamics simulations could be explained by the fundamental-adjoint pure gauge 
transition, and the real cause of this phenomenon remains to be uncovered. 
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